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Abstract 

Let L be the EucHdean functional with p-th power- weighted edges. Examples 
include the sum of the p-th power- weighted lengths of the edges in minimal 
spanning trees, traveling salesman tours, and minimal matchings. Motivated 
by the works of Steele, Redmond and Yukich (1994, 1996) have shown that 
for n i.i.d. sample points from [0,1]^, . . . , X„})/n(^-P)/'^ 

converges a.s. to a finite constant. Here we bound the rate of convergence of 
SL({Xi,...,X„})/n('^-P)/'^. 
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1 Introduction. 

Let {Xi, . . . be n i.i.d. sample points from R'^, d > 2, and let < p < oo. A 
traveling salesman problem (TSP) is to find a permutation vr on {1, . . . , n} such that 



(i+i) ~ ^Hj) 



IP 



= min < — A^/(j)|^ : n' a permutation on {1, ... , n} 

where |Xj— Xj| is the Euclidean distance between Xj and Xj and where 7r(n-|-l) := 7r(l) 
and Tr'^n + 1) := 7r'(l). Let Lt5p ({Ai, . . . , A„},p) be the sum of the p-th power- 
weighted lengths of the edges in a minimal tour vr. In the case {Ai, . . . , A„} = define 
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Ltsp{^,p) = 0. Beardwood, Halton, and Hammersley (1959) showed that there exists 
a strictly positive but finite constant a{LTSP, d, 1) such that for i.i.d. sample points 
{Xi : i > 1} with common distribution /i, which has a compact support in R*^, d > 2, 
as n — > oo 



a{LTSP, d, 1) J f^''-^^'\x)dx a.s. (1.1) 



where / is the density function of the absolutely continuous part of /x. 

The asymptotic behavior (1.1) of the TSP functional is not an isolated one. A 
minimal matching (MM) on {Xi, . . . , X„} is a permutation tt on {1, . . . , n} such that 



'-7r(2j-l)| 



[n/2] 




E 










f [«/2] 


min 









where [n/2] is the largest integer I with I < n/2. Let Lmm ■ ■ ■ ,X„},p) be the 
sum of the p-th power-weighted lengths of the edges in a minimal matching tt. In the 
case {Xi, Xn} = define Lmm{^,p) = 0. 

A minimal spanning tree (MST) on {Xi, . . . , Xn} is a spanning tree T on the given 
point set {Xi, . . . , X^} such that 



{Xi,Xj)eT 



\p 



— min < ^ \Xi — Xj\^ : T' a spanning tree on {Xi, . . . , Xn} > . 
[iXi,Xj)eT' J 

Let Lmst {{Xi, . . . , X„},p) be the sum of the p-th power-weighted lengths of the edges 
in a minimal spanning tree T. In the case {Xi, . . . , Xn} — define Lmst{^,p) — 0. 

A Steiner minimal spanning tree (SMST) on {Xi, . . . ,Xn} is a spanning tree T 
on a point set containing (we call such T a Steiner spanning tree on 

{Xi, . . . , Xn}) such that 



{Xi,Xj)eT 



.IP 

3\ 



— min < ~ Xj\^ ■ ^' ^ Steiner spanning tree on {Xi, . . . , Xn} > . 

{Xi,Xj)eT' 
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Let Lsmst{{^i, ■ ■ ■ ,^n},p) be the sum of the p-th power- weighted lengths of the 
edges in a minimal Steiner spanning tree T. In the case {Xi,...,Xn} — define 
LsMsri^-iP) — 0. Note that for p > 1, Lsmst — 0. So, whenever we talk about Lsmst, 
we always consider the case < p < 1. 

A rectihnear Steiner minimal spanning tree (RSMST) on {Xi, . . . , Xn} is a Steiner 
spanning tree T on {Xi, . . . , Xn} in which all the edges are rectilinear (we call such T 
a rectilinear Steiner spanning tree on {Xi, . . . , X^}) such that 

{Xi,Xj)eT 

— min < \Xi — Xj\'^ : T' a rectilinear spanning tree on {Xi, . . . , Xn} > . 

(iXi,Xj)eT' J 

Let LjiSMST {{^1: ■ ■ ■ :Xn}:P) be the sum of the p-th power-weighted lengths of the 
edges in a minimal rectilinear Steiner spanning tree T. In the case {Xi, . . . , Xn} — 
define Lrsmst{^-iP) — 0. Note again that for p > 1, Lrsmst — 0. So, whenever we 
talk about L^smst, we always consider the case < p < 1. 

In a series of papers Steele (1981A, 1981B, 1988, 1990) showed that the asymptotic 
behavior (1.1) appears for various functionals including some of the above five. Red- 
mond and Yukich (1994, 1996) and Lee (1999) further developed the general conditions 
providing the asymptotic behavior (1.1). One may consult Section 1.2 of Yukich (1998) 
for a brief history of this field. 

Our results are stated below. But, first we would like to spell out the restrictions 
on the Euclidean functional L. Wc call L{A,B,p), A a finite subset of a box B = 
nf=i[^«'^« + '^]' ^ ~ (^1) • • • j^d) £ IR^, < s < oo, d > 2, < p < oo, a subadditive 
Euclidean functional (or a weak Euclidean functional) of power p if the following four 
conditions are met: 

L(0,5,p) = O, (1.2) 

for y eW'' and < t < oo 

L{y + tA,y + tB,p)=t^L{A,B,p), (1.3) 

\L{A, B,p) - L{B, B,p)\< C\AAB\^'^-P^/'^sP, (1.4) 
and for a partition {Qi,l < i < m*^} of [0, 1]*^ into m'' subboxes of edge length 

L(^,[o,i]'^,p) <E^(-4ng,,gi,p)+Cm'^-^. (1.5) 



Rates of convergence of means of Euclidean functionals 



4 



By (1.2) and (1.4) with s = 1 and B — we have for a finite subset A of the unit 
box [0, 1]'^ and for d > 2, < p < d, 

|L(^,[0,l]^p) <C|>l|('^-f)/<^. (1.6) 

More strongly, for all the Euclidean functionals L of our interest in this paper (look at 
Theorem 2 below for the full list of such L) based on the space fiUing curve heuristic, 
as shown by Steele (1997), there is an extension of the above bound (1.6) which covers 
both the d — 1 case and the p > d case; for a finite subset A of the unit box [0, l]'' and 
for d > 1, p > 0, 

\L{A, [0, lY,p)\ < C{\A\^''-P^/'' V 1). (1.7) 

Note that (1.6) follows from the assumptions of the Euchdean functional whereas (1.7) 
follows from the specific feature of the Euclidean functional of our interest in this paper. 
Let Un be n i.i.d. uniform points in [0, l]**. 



Theorem A. [Redmond and Yukich (1994, 1996), Lee(1999)] Let L be a weak Eu- 
clidean functional of power < p < d. Then there exists a finite constant a :— 
a{L, d, p) such that as n ^ oo 

L(Wn,[0,l]^p) 



a c.c. and in L^, (1.8) 



where ^ a c.c. (complete convergence) means that for any £ > 0, P(|F„ — a| > 

e) < oo. 

For a typical weak Euclidean functional L, the limit a in (1.8) is strictly positive: 
In most situations of interest, the limit a is just the subadditive constant and therefore 
must be strictly positive. 

We call L*{A, B,p) a superadditive Euclidean functional of power p if —L*{A, B,p) 
is a subadditive Euclidean functional of power p. We call L a Euclidean functional of 
power p if L is a subadditive Euclidean functional of power p, if L* is a superadditive 
Euclidean functional of power p, and if 

L*{A,[0,ir,p)<L{A,[0,lf,p), (1.9) 

and for the n uniform points W„ in [0, 1]*^ 



\EL{Un, [0, l]^p) - EL*{Un, [0, l]^p)| = oin^'-P^/"). 



(1.10) 
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Theorem B. [Redmond and Yukich (1994, 1996), Lee(1999)] Let L be a Euclidean 
functional of power < p < d. Then for i.i.d. sample points {Xi : i > 1} with common 
distribution which has a compact support in [0, l]'', as n — > oo 

L({X„.. ,X } [0,inp) r ^^ 

where a :— a{L,d,p) is a finite constant given by (1.8) and where / is the density 
function of the absolutely continuous part of /i. 



If there is no confusion, to save the heavy notations from now on in the case B = 
[0, l]'^ we use the notation L{A) and L*{A) instead oi L{A, [0, 1]"^,^) and L*{A, [0, 1]"^,^), 
respectively. 

We call L a strong Euclidean functional of power 0<p<(iifLisa Euclidean 
functional of power p and if for the homogeneous Poisson point process Vn of intensity 
n on [0, 1]'^ 

^^{d-i-p)/d for < p < - 1 
C ior p = d — 1 = 1, 



Cn^d-i-p)/d for < p < d - 1 

C(lognVl) forp = ci-l7^1 

C ior d-1 <p < d 

C for p = d — 1 = 1, 



|£;L(ZY„±fe)-£;L(w„)| 

( C'(A;('«-i-p)/('«-i) VI) for 1 < A; < n^'^"^)/'^ 

< \ Ckn^P/'^ for n^^-^^^'^ < k < n/2,0< p < d ~ 1 (1.14) 
( Ckn-^'^-^^/'^ for n^'^-^^/'^ <k< n/2, d - 1 < p < d. 

We call L a very strong Euclidean functional of power 0<p<difLisa strong 
Euclidean functional of power p and if L satisfies the add-one bound; 

\EL{Un+i) - EL{Un)\ < Cn-P'''. (1.15) 



Our first result is as follows. 
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Theorem 1. (i) For a strong Euclidean functional L of power < p < d, 



Q^{d-l-p)/2{d-l) 



(j^{d-i-p)/d for < p < - 1 

C(lognVl) forp = d-l7^1 

C ioY d-l <p <d 

C for p = d - 1 = 1, 



-C 
-C 
-C 



< EL{Un) - an' 



{d-p)/d < 



where a := a{L,d,p) is a finite constant given by (1.8). 
(ii) For a very strong Euclidean functional L of power 



< p < d, 



Cn^--d < EL{Un) - an^'^-P^''^ < 



I 



C(logn V 1) 
C 

c 



for < p < (i - 1 
ioi p = d — 1 ^ 1 
ioi d — 1 < p < d 
lox p — d — 1 — 1, 



where a := a{L, d,p) is a finite constant given by (1.8). 

To see the point of Theorem 1 we compare it with Theorem 5.2 of Yukich (1998). 
In Theorem 5.2 of Yukich (1998), he showed that, if L satisfies the close in mean 
approximation, that is, \EL{Un) — EL*{Un)\ = o(n('^~^)/'^), and the add-one bound, 
then 



The provided rate is quite satisfactory. However, as we see in Theorem 2 below, many 
typical L do not satisfy the add-one bound and hence we cannot apply Theorem 5.2 of 
Yukich (1998) to those L. We wish to provide the same rate without the add-one bound 
so that typical L has the provided rate of convergence. What we find in Theorem 1 is 
that this task can be done by strengthening the close in mean approximation to (1.12) 
and by adding two extra conditions (1.13) and (1.14). We also find that for the case 
p — d — 1 ^ 1 typical L do not satisfy the other condition used in Theorem 5.2 of 
Yukich (1998) and we properly fix it in (1.12) and (1.13). We also take care the case 
< p < 1 which was excluded in Theorem 5.2 of Yukich (1998). 

Note that the bounds in Theorem 1 follow from the assumptions of the Euclidean 
functionals. If we use the specific property of the Euclidean functional, in some cases 
we can get much better bounds and surprisingly we can even get a strictly positive 
lower bound for large n; see Jaillet (1993), Rhee (1994), Yukich (1998), Lee (2000) for 
these specific results. 

For a given subadditive Euclidean functional L, to find the rate of convergence of 
ELiUn) using Theorem 1 we have to construct a superadditive Euclidean functional 
L* which puts us into the situations considered in Theorem 1. There may be several 



EL{Uu) — an' 
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ways to construct such an L*. The successful L* is the one which uses the boundary 
freely. For example let's consider the traveling salesman problem. Suppose that there 
are n cities A in [0, l]''. Then Ltsp{A, [0, l]'', 1) is the total mileage to travel all the n 
cities. Suppose however that there are "free" ways along the boundary d[0, 1]^ of [0, l]'' 
in which the government pays the gas. In this case we may save some gas by traveling 
along the boundary and L^sp(A, [0, l]'', 1) is the total mileage. 

In the case 1 < p < d following the idea of Redmond and Yukich (1994, 1996) 
we construct superadditive Euclidean functionals L* for the MM, MST, TSP, SMST, 
RSMST (of course for the SMST, RSMST we consider the p = 1 case only). They are 

L*mm{A: B,p) :— min {Lmm{A \JB,B,p) : B a, finite subset of dB} , 
L*mst{A,B,p) Lmst{A,B,p) Amm!^y^^LMST{Aj[J{bj},B,p)^ , 
where the minimum is taken over the partition {.4^} of A and bj e dB, 

L*j,sp{A,B,p) := Ltsp{A,B,p) A min I^Y^LrspiAj U {6,-, fc^}, S,p) | , 

where the minimum is taken over the partition {Aj} of A and bj, b'j e dB and where 
for a finite subset {Xi, X^} of B with \{Xi, X„} ndB\>2 



' n—l 



Ltsp{{Xi, . . ■,Xn}, B,p) := min |^ |X7r(j+i) - J , 

where the minimum is taken over the permutation tt on {1, . . . ,n} such that ^7r(i)) 
^T^in) £ dB. Note that in the definition of Ltsp the sum is up to n — 1 so that we 
travel free from X7r(„) to X7r(i) along the boundary dB. Similarly, 

L*smst{A,B,p) := Lsmst{A,B,p) A min |^-L5MST(^j U {6j},B,p)| , 

where the minimum is taken over the partition {.4^} of A and bj e dB, 

L*rsmst{A, B,p) := Lrsmst{A, B,p) A min | J]] -Lr5mst(A- ^ i^j)^ B,p)^, 

where the minimum is taken over the partition {Aj} of A and hj G dB. In the case 
^ = define L^^(0,S,p) = L;,5^(0,S,p) = L*.,sp{%, B ,p) = LJ^5^(0,5,p) = 
L*rsmst{^^B,p) = Q. 
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In the case 1 < p < d following the idea of Lee (1999) we construct superadditive 
Euchdean functionals L* for the MM, MST, TSP, SMST, and RSMST. They are given 
in the following way. In the above L* for 1 < p < d, there are some edges {X, Y) from 
a boundary point X e dB. In the case < p < 1, for any matching, tree, or tour, we 
don't pay the full price for the edge {X, Y) from the boundary point X e dB; for this 
edge we pay half of the full price \X — y|^, i.e., \X — Y\^ /2. 

Theorem 2. (i) For the Euchdean functional L of MM, MST, TSP with <p < d 
and SMST, RSMST with < p < 1, L is a strong Euclidean functional (and hence 
Theorems 1 (i), 3, 4 can be applicable to this L). 

(ii) For the Euclidean functional L of MST with < p < and TSP, SMST, 
RSMST with 0<p< l,Lisa very strong Euclidean functional (and hence Theorems 

1 (ii), 3, 4 can be applicable to this L). 

In Section 2, we develop a theory on the rate of convergence of EL for the uniform 
sample points, i.e., we prove Theorems 1 and 2. In Section 3, we continue to build a 
theory on the rate of convergence of EL for the non-uniform sample points which was 
started by Hero, Costa, and Ma (2003). 

2 Rates of convergence; the uniform case 

In this section we prove Theorems 1 and 2. The main idea comes from the symmetry 
argument of the patching in Alexander (1994). Using this symmetry argument we get 
the nice moment estimate. 

Lemma 1. Let L be a strong Euclidean functional of power < p < d. Then, 




for < p < d - 1 

for p = d — 1 ^ 1 
for d — 1 < p < d 
for p = d — 1 = 1. 



where a := a{L,d,p) is a finite constant given by (1.8). 
Proof. By the usual subadditive argument (see page 54 



of [16]) for L, 



C < EL{Vn) - an' 



{d-p)/d 



By the same argument for L*, 



^L*(P„) -an('^-f)/'^ < C. 
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Now, use (1.13) to get the Lemma. | 

Lemma 2. Let Nn is a Poisson random variable with mean n which is independent of 
i.i.d. uniform points {Xf, 1 <i < oo} on [0, 1]'*. 

(i) If L is a strong Euchdean functional of power < p < d, then 

\EL{Un) - EL{Un^)\ < v 1). (2.1) 

(ii) If L is a very strong Euchdean functional of power < p < then 

\EL{Un) - EL{UnJ\ < Cn^/2-p/d_ ^2.2) 

Proof, (i) If L is a strong Euchdean functional of power < p < d — 1, by (1.14), 
(1.6) and by Jensen's inequality we have 

\EL{Un) - EL{UnJ\ < \E{L{Un) - L{UnJ)1 {\Nr,-n\<n(d-i)/dy\ 

+ \E{L{Un) - L{Un„))1 

{iV„<n/2}| 
{iV„>3n/2} I 

< CE{\Nr, - V 1) + Cn-Pl'^E\Nr, - n\ 
+Cn^^-p)/dp(N^ < n/2) + CEN^^-^^'%^^^.,,,,,} 

< C(n(''-^-f)/2('^-^) V 1) + 

+Cn^'i-p)/<^P^N^ < n/2) +CEiV('^-P)/'^l{^„>3n/2}. 

Since Nn has a very light tail, the last two terms of the left hand side in the above 
inequality are negligible and (2.1) follows. The argument for d — 1 < p < dis same. 

(ii) If L is a very strong Euclidean functional of power < p < d, there is a standard 
argument from Theorem 5.2 of Yukich (1998): By (1.6) we have 

\EL{Un)-EL{UNj\ < Cn-P'''E\Nn — l{l<|Ar„-n|<n/2} 

+^^(rf-p)/rfp(^^ < n/2) +CEiV('^-^')/'^l{^„>3n/2}. 

Thus, (2.2) follows. | 

Proof of Theorem 1. Theorem 1 follows from Lemmas 1 and 2. | 

Lemma 3. For the superadditive Euclidean functional L* of MM, MST, TSP with 
< p < d and SMST, RSMST with < p < 1, let NsiJ^n) be the number of points in 
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Vn which are connected to the boundary (9[0, l]*^ in the minimal graph which is used 
to calculate L*. Using the same minimal graph let LsiVn) be the sum of the p-th 
power-weighted lengths of the edges connecting points in Vn to the boundary 9[0, l]''. 
Then, 

ENsiVn) < Cn('^-^)/^ 
(j^(d-i-p)/d for < p < d - 1 
C(lognVl) forp^d-l^^l 
C ioT d-1 < p < d 

C for p = (i — 1 = 1. 

The same estimates for W„ instead of Vn also holds. 

Proof. We just follow the argument for Lemma 3.8 of Yukich (1998) or Lemma 4 of 
Lee (1999) and dig out the quantities of interest. We skip its proof. | 

Lemma 4. For the subadditive Euclidean functional L of MM, MST, TSP with 
< p < d and SMST, RSMST with < p < 1, 

for < p < d - 1 

0<EL(P„)-EL-(P„)<^ 'ZTJ~<ltl P.3) 

C for p = d - 1 = 1. 

The same estimates for Un instead of Vn also holds. 

Proof. Since the case 1 < p < d has been known quite sometime (see Lemma 3.10 of 
Yukich (1998)), we consider the case < p < 1 only. Even in this case, the argument 
for the case l<p<d— lof Lemma 3.10 of Yukich (1998) works; we just need to use 
the corresponding estimates for < p < 1 in Lemma 3. | 

Lemma 5. For the subadditive Euclidean functional L of MM, MST, TSP with 
< p < d and SMST, RSMST with < p < 1, L satisfies (1.14). 
Proof. Since all the arguments are similar, here we prove the Lemma for the MST 
case only. We leave all the other cases to the reader as an exercise. First we claim that 
for 1 < A; < 

EL{Un+k) < ELiUn) + C(A;('='-i-^')/('^-i) V 1). (2.4) 

Let {Xi, . . . ,Xn+k} be the n + k i.i.d. uniform points on [0, l]'^. By renaming these 
points Xi := (X/, . . . , Xf) we may assume that the first coordinates of these n+k points 
are increasing, i.e., Xl < X^ < ■ ■ ■ < X^+j^.. With this renaming, let t :— 1 — X^+i. 
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Figure 1: A way to construct a spanning tree on {Xn+i, ■ ■ ■ , Xn+k}', add an edge 
{Xi, Xj) if an edge {Yi, Yj) is in tlie MST on {Fn+i, • • • , ^n+fe}- 



After constructing an MST on {Xi, . . . , Xn} and another MST on {Xn+i, . . . , 
by adding an edge {Xn, Xn+i) we have a spanning tree on Un+k- So, we have 

EL{lln+k) < EL{{Xi, . . . , Xn}) + EL{{Xn+l, . . ., Xn+k}) + C. 

Let X* := ((1 - r)-^Xj,X], Xf, . . . , Xf),l < j < n. Then, since {X^, ...,X*} are n 
i.i.d. uniform samples from [0, 1]'^, 

EL{{X,, ...,Xn})< EL{{X*, ...,X*n})^ EL{Un), 

and hence 

EL{Un+k) < EL(Un) + . . . , Xn+k}) + C. (2.5) 

Now, we estimate EL{{Xn+i., ■ ■ ■ , Xn+k})- Let 1^, n + 1 < i < n -\- k, he the 
projection of X^ to {1} x [0, 1]^~^. We construct an MST on {Yn+i, . . . , Yn+k} and then 
using this MST we construct a spanning tree on {Xn+i, ■ ■ ■ , Xn+k} by adding an edge 
{Xi, Xj) if an edge {Yi, Yj) is in this MST (see Figure 1). For this edge {Xi, Xj) 

- Xj\^ < {\Yi - Yj\ + 2Tf < C\Yi - Yj\P + CtP. 
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Figure 2: A way to construct a connected graph on Bf, add an edge (X^.X^/), Xi G Bi, 
Xi' G Bi, if {Xi,Xj), Xj G Wj, and {Xi',Xj'), X-i G >Vi, were edges in the original 
MST T{Un) and if {Xj, ) is an edge in the MST T{yVi). 

Since £;t*' < C{k/n)P, we have then by (1.7) 

EL{{Xr,+,,---,Xn+k}) < CEL{{Y^+,,...,Y^+k},[0,n''~'x{n,P)+CkETP 

< C'(A;(<i-i-f)/(<i-i) V 1) + CkErP 

< c(k^'^-^-p^/^'^-^^ y 1). (2.6) 

Therefore, (2.4) follows from (2.5)- (2.6). 

Second, for n^'^~^y'^ < k < n/2 we iterate the above argument [k/n^'^~^^^'^] + 1 times 
and we have 

EL{Un+k) < EL(Un)+C([k/n^''-^^/''] + l)(n^'-'-P^/''yi) 

< FT(ir\^l forO<p<d-l 

< EL{Un) + I ckn-^^-')l^ for d - 1 < p < d. ^^"^^ 

Third, we claim that for 1 <k < n^'^~^'>/'^ 

EL{Un-k) < EL{Un) + C{k^'^-^-P^'^'^-^^ V 1). (2.8) 

As we did before, by renaming the n random points lAn '■= {Xi, . . . , X„} we may assume 
that the first coordinates of these n points are increasing, i.e., X^ < X^ < ■ ■ ■ < X^. 



Rates of convergence of means of Euclidean functionals 



13 



For 1 < i < n - k + 1, let Bi = Un \ {Xi, Xi+k_i} and let n := X^^f, - X/. We 
construct an MST T{Un) on Un- Prom the MST T{Un), we remove all the points Xj, 
i<j<i + k — 1, and all the edges in the MST T{Un) which use those Xj as one 
end. The resulting graph Di on Bi is disconnected. By adding some extra edges to 
Di we will construct a connected graph on B^. To do this we first collect those Xj, 
i < j < i + k — 1, ior which Xj was connected to a remainder point Xk G Bi by 
an edge in the original MST T{Un)- We call the set of all those collected Xj as Wj. 
We construct an MST T{Wi) on Wi. Using this MST T{Wi) we have the following 
addition rule of edges to the disconnected graph on Bi] we add an edge {Xi,X,^'), 
Xi e Bi, X^ G if Xj e Wi, and (X;.,Xy), X.> e Wi, were edges in the 

original MST T(W„) and if {Xj,X.') is an edge in the MST r(Wi). By adding these 
extra edges to Di we get a connected graph on Bi (see Figure 2). 
In this case 

\Xi - Xi>\P < C\Xi - Xj\P + C\Xj - Xj,\P + C\Xj, - Xi,\p. 

Since the degree of a vertex in an MST on the set of points in M.'^ is bounded by a 
universal constant which depends only on the dimension d, by the argument of (2.6) 

EL{Bi) 

i+k-l 

< EL{Un) + C J] E\Xj - + CEL{Wi) 

j=i {Xj,Xk)eT(Un) 

i+k-l 

< ELiUn) + C Y - ^'^1'' + C{k^'^-^-^^'^'^-^^ V 1). (2.9) 

j=i (X,-,Xfc)GT(W„) 

Let X* := Xj — TjCi for i + k < j < n, where ei := (1,0,..., 0)) is the first unit vector, 
and let B* := {Xi, . . . , Xj_i, X*^^, . . . , X*}. B* then consists of the n-k i.i.d. uniform 
samples from [0, 1]'^"^ x [0, 1 — Tj]. By its construction, 

LiB*) < L{Bi). 

Multiplying the first coordinates of the points of B* by (1 — Ti)~^ we have n i.i.d. 
uniform samples from [0, 1]"^. Hence, since (1 — Tj)~*' < 1 + Ctj for the case Tj < 1/2, 
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and since 1 + Ctj is increasing in Tj and E{L{B*)\Ti) is decreasing in Tj, 

EL{l/{n-k) 

< E{{1- n)-PE{L{B*)\n)\n < 1/2) P{n < 1/2) + E {L{Un-k)-,n > 1/2) 

< E{{i + Cn)E{L{!Sl)\n)\n < 1/2) p(t, < 1/2) + 6-^^^ 

< {E{l + Cn\n < 1/2) E{E{L{B*)\Ti)\Ti < 1/2)) P(t, < 1/2) + e"^'^ 

< (l + E {E{L{B;)\n)\n < 1/2) + e-^" 
k \ EL{Bt) , 



< {i + c-\EL{Un) + c E\x^-Xk\'' + c{k^'^-^-^^'^'^-^^yi). 

^ (x,-,Xfc)eT(w„) 
Since k < n'-'^"^^/'^, we have 

-EL(Ur,) < C-n^"-^^^" < Ck^d-i-pVid-i) _ 
n n 

So, by combining the terms C^EL{Un) and C{U<^-'^-p^/^'^-'^'^ V 1) and by increasing the 
constant C in the term C(A;^'^-^-^')/('^-^) V 1) we have 

i+k-l 

EL(Un-k) <EL(Un)+C J2 Yl E\Xj-Xk\P + C(k^'^-^-P^/^'^-^^yi). (2.10) 

j=i {Xj,Xk)eT(Un) 

Here comes the highhght of the argument; we use the symmetry argument to estimate 
the term E^"^ E(x„x,)eT{w„) ^ - ^k\^- Since 

n—k+l i+k-l 

J2 Yl E E\Xj-Xkf <CkEL{Un), 
i=i j=i {Xj,Xk)eT(Un) 

averaging (2.10) over l<i<n — k + 1 (by (1.6)) we have (2.8). 

Last, for n^^~^^^'^ < k < n/2 we iterate the above argument \k/nS'^~^'>/'^'^ + 1 times 
and we have 

EL{Un-k) < £;L(W„) + C([A;/n('^-^)/'^] + l)(n('^-^-^')/'^Vl) 

< EL{Un) + I for d - 1 < p < d. (^--^^^ 



Rates of convergence of means of Euclidean functionals 



15 



I 

For the Euclidean functional L of MST with < p < d, L satisfies the add-one 
bound (1.15) as shown by Redmond and Yukich (1994). In Lemma 6 below we show 
that for the case < p < 1 many typical L also satisfy the add-one bound (1.15) 
and hence we provide some affirmative answers to the issue raised in p. 55 of Yukich 
(1998). However, we cannot handle the case 1 < p < d for those L and more seriously 
we cannot prove the add-one bound for the minimal matching Euclidean functional 
Lmm- So, we think the add-one bound (1-15) condition is very restrictive. 

Lemma 6. If L is either the MST Euclidean functional L with < p < or the TSP, 
SMST, RSMST Euclidean functional with < p < 1, then L satisfies (1.15). 
Proof. Since all the arguments are similar, here we prove the Lemma for the TSP case 
only. We leave all the other cases to the reader as an exercise. 
Fix < p < 1. First we claim that 

EL{Un) < EL(Un+i) + Cn-'^l'^. (2.12) 

There are two edges (X^, Xj) and (Xj, X^) adjacent to Xj in the minimal tour T(W„+i). 
Remove these two and add the edge (Xj,Xfe). Then, we have a tour on \ {Xj}. 
Since 

~ ^ ~ ~ -'^fel)^ ^ C |Xj — Xj|^ -|- C |Xj — Xfcl^ , 

we have 

EL(W„+i \ {X,}) < EL{Un^^) + C ^ - ^il'' • 

Since EL{pin+i \ {X^}) = EL{Un), and since El'i^ E(Xi,x,)er(w„+i) ^ 1^* - ^jf = 
2EL{Un+i) < Cn^'^~P^^^ (by (1-6)), averaging the above inequality over 1 < i < n + 1 
we have (2.12). 

Now, we claim that 

EL{Un+i) < EL{Un) + Cn-Pl'^. (2.13) 

Let T(Un) be the minimal tour on W„. For a given X„_|_i, among n i.i.d. uniform points 
Un find a nearest point Xj to X„+i and let (Xj.X,) be an edge in the minimal tour 
T(Un)- By removing this edge and by adding (Xj,X„_|_i) and (X„+i,Xj) we construct 
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a tour on Un+i- Since < p < 1, we have < B^ + C' for a triangle with side lengths 
A, B, C. So, using this triangular inequality 

L{Un+i) < — \Xi — Xj\^ + \Xi — + \Xj — 

< L{Un) + 2 \Xi — Xn-\-lf 

= L{Un) + 2 min \Xk - X^+^l^ 

l<k<n 

By taking expectations in the above inequality we have (2.13). | 

Proof of Theorem 2. (i) follows from Lemmas 4 and 5 and (ii) follows from Lemmas 
4, 5, and 6. | 



3 Rates of convergence; the non-uniform case 

In this section, we continue to build a theory on the rate of convergence of EL for 
the non-uniform sample points which was started by Hero, Costa, and Ma (2003). In 
Hero, Costa, and Ma (2003) they studied the rate of convergence of EL for the non- 
uniform sample points based on the very restrictive add-one bound (1-15) condition. 
Our starting point of this study is to try to remove this condition in their argument 
and our work in Section 2 is the result of this trial. In this section we build a theory 
on the rate of convergence of EL for the non-uniform sample points which does not 
depend on the add-one bound (1-15) condition. 

First, we work with a block density function; a probability density function cf) of 
the form 

i=l 

where 0^ > is a constant and where {Qi, I <i < m'^} is a partition of [0, 1]*^ into 
subboxes of edge length ■m~^, is a block density function of level m. 

Theorem 3. Let L be a strong Euchdean functional of power < p < d. Then for 
i.i.d. sample points {Xi . i > 1} with common block density function of level m there 
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exists a constant C > 0, independent of (j), n, m, such that 
EL{{X^,...,Xr,}) 



< 



C{nm-'^)-^/'^ for < p < (i - 1 

C(nm-^)-('^-P)/'^ (log(nm-'^) V l) for p = d - 1 ^ 1 , . 

C(nm-'^)-('^-f)/'=' for d - 1< p < d ^ ^ ^ 

clnm-'^y^'^-py'^ for p = d - 1 = 1, 

where a is a constant given by Theorem A. 

Proof. Hero, Costa, and Ma (2003) consider the case 1 < p < d — 1 only. However, 
their argument also works for all the other three cases considered in this theorem 
without any major changes. For reader's convenience here we reproduce the argument 
of Hero, Costa, and Ma (2003) for the case d — 1 < p < d and the case p = d — 1 7^ 1 
with appropriate (minor) changes. 

First, we handle the case d—1 < p < d. Let Xi, . . . , X„ be n i.i.d. samples with the 
common block density function of level m and let L be a strong Euclidean functional. 
Let rii be the number of samples {Xi, . . . falling into Qi. Then, by (1.5), (1.3), 
Theorem 1 (i) 

EL({Xi,...,Xj,[0,l]^p) < 5^EL({Xi,...,X„}nQ,,g,,p) + Cm'^-f 

i=l 

= m-Pj2EE(L({Ui, Un,}, [0, Ij^p)]^) + Cm^-P 

i=l 
i=l 



i=l 

d 



By Jensen's inequality, for < 1/ < 1 and pi := cpim 



So, 

ELi{X,,...,X^},[0,ir,p) .U^-Pyd^-d , c 



^{d-p)/d — * {nm~'^Y'^~py'^ 

(j)'^^-py\x) dx + 



1=1 

C 



— a 



[0,1] 
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Similarly, by (1.12) 

EL*{{X^,...,Xn},[0, l]',p) 
> m-Pj^EE {L*{{U,, UnJ, [0, lY,p)\ni) - Cm^-^ 



i=i 



> m-P J2 EE (l({C/i, . . . , C/„J, [0, l]^p) - C|n,) - Cm"- 

i=l 

> m-Pj^E (anf -^^/'^ - C - c) - Cm^'P 



i=l 



rij \ (.d-p)/d 



Now, we claim that for < < 1 and Pi := (pifn " 



E[^] >pt-p, 



-1/2 



(3.2) 



If g e C^(0, oo) and if g is concave over x > 0, monotone increasing over x > 0, and 
g'(O) = 0, then for any xq > 0, 

g{x) > g{xo) - - Xo\. 

Xq 

Thus, with g{x) := x^, < u < 1, x := Ui/n, and xq := pi, we have 



Hi 



n 



Tii 



n 



Pi 



Take the expectation on both sides. Since by Chebyshev's inequahty 



E 



Pi 

n 



< E 



2\ 1/2 



< 



indeed wc have (3.2). 
So, by (3.2) 



EL*{{X,,...,Xn},[0,l]'',p) 



> a [ ct>^'^-P^/\x)dx ^ / 05-§(x)dx- 



C 
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Therefore, by (1.9) and (1.12) we have 

EL{Xi, . . . , Xn) 



< 



n 



a 



(d-p)/d 



— a 



'[0,1]' 



02 d(^x)dx + 



C 



Second, we handle the case p = d— 1^1 in a. similar manner. 
EL{{X,,...,X^},[0,lY,p) 

< m-Pj2^ (anf "^^/'^ + C (log V 1)) + Cm'^-P 



i=l 



< am-Pn^'^-P^/'^ ^ E (^— ) ^'^ ^^^"^ + Cm-P ^ E (log rii) + Cm'^-P . 

i=l ^ i=l 

Since rii is highly concentrated around its mean npi, for a large but fixed C we have 
E^lognj < Clognpi. So, 

£;l({Xi,...,x4,[o, i]^p) 



< a 



E 

i=l 



dd-p)/d -d . <^E£'i(log0»)"^~'^ , Clog(nm-'^) C 



7 



C L^.Jog(^{x)dx c\o 



(nm~'^) C 



^ °^ l[o,i]'i^ (a;)(ia;+ ^^^_dyd-p)/d ^^^-dyd-p)/d ^ ^^^-dyd-p)/d- 



By the same reasoning, 

£;l-({Xi,...,x4,[o, i]^p) 

7[o,i]d (nm 
Clog(nm-'^) C 



> a 



a f i_p <^/[o,i]rflQg0(^)^^ 

(^^-<i)(<i-p)/d 



Therefore, by (1.9) and (1.12) we have 
EL{Xi, . . . , X„) 



< 



^{d-p)/d 

c 



— a 



-p)/d, 



x) dx 



[0,1]" 



(nm-'^)(^-p)/<^ 



/ log(j){x)dx+ ^TjT^ / 



)2 d(^x)dx 



C 



Clog^ 
{nm~'^) i'i-p)/d [nm-<^) id-p)/d ' 
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I 

Second, we work with a probability density function / in the Holder class '^{P, K, [0, 1]*^), 
/3 > 0; a probability density function / : [0, l]*^ ^ R is in K, [0, 1]'*) if for 

x,y & [0, 1]*^ 

\m-pi^^{y)\<K\x-yf, 

where |_/3J is the largest integer I with I < /3 and where plf^ is the Taylor expansion of 
/ at X up to degree [P\ . 



Lemma 7. If a probability density function / is in [0> M^)^ < /3 < 1, then 

for any partition {Qi, 1 < i < m'^} of [0, 1]'' we define an approximating block density 
function with level m of the form = YliLi by (pi — fn'^ Jq. f{x)dx. For 

this approximating block density function 0, 



/ 



|/(a;) - (j){x)\ dx < df^'^Km-^. (3.3) 

'[0,1]'' 

Proof. By the mean value theorem there exists a point r]i e Qi such that 

(t>i^rn'^ / f{,x)dx = f{r]i). 
JQi 

With this rii, 

/ \4>{x) - f{x)\ dx ^ J2 - 

i=i JQi 



IIL n 

< / K \rii — xf dx 

i=i JQi 



< d^/''Km-f. 



Theorem 4. Let L be a strong Euclidean functional of power < p < d. Then 
for i.i.d. sample points {Xj : i > 1} with common probability density function / e 
Xl(/9, K, [0, 1]'^), < /3 < 1, there exists a constant C := C{d, P,p) > such that 



/ fj(d~p)/d 
Cn {p{d-p)/d+i)d 



EL{Xi, . . . , Xn) 



^{d-p)/d 



a 



dx 



< I 



Cn (logn V 1) 

Cn o+d 

I3(d-p)/d 

Cn i^+d- 



for < p < d - 1 
for p = d — 1 ^ 1 
ior d — 1 < p < d 
ior p — d — 1 — 1, 
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where a is a constant given by Theorem A. 

Proof. As we did in Lemma 7, for any partition {Qi, 1 < i < m^} of [0, l]"^ we define 
an approximating block density (j) with level m of the form = YliLi4>i'^Qii^) by 
0j = m'^ Jq f{x) dx for which we have (3.3). 

Here we consider the case < p < d — 1 only. We leave all the other cases to 
the reader as an exercise. Let {-^^1}^=! be i.i.d. with the common probability density 
function (f). Then, 

EL{{X,,...,Xn}) 



< 



EL{{X,,...,X^}) 



-a I f^'^-P^'\x)dx 



a / (j)^'^-P^I'^{x)dx 



+Q! 



[0,1] d 



{d-p)/d. 



x)dx — 



-P^/'^ix) dx 



[0,1]" 



+ 



i?L({Xi, . . . , XJ) - EL({Xi, . . . , XJ) 



:= 7 + 77 + 777. 



j^{d-p)/d 

The term I can be handled by Theorem 3 and it is bounded by C{nm~^)~^^'^. 
By (3.3), the term II is bounded by Cm~^^'^~P^^^; 



II < a 



< a [ \f{x) - 0(3 

[ \f{x)-cl>{x)\dx' 



yx^^"-"^'" dx 



< a 



(d-p)/d 



^'^-P^/^ ^-f3(d-p)/d 



By standard coupling and by (1.6), the term III is bounded by CmT^^'^'P^^'^] 



III = 

< CE 



EL{{Xi, . . . , XJ) - EL{{X,, Xn}) 

(d-p)/d 



{Xi, . . . , X„} A{Xi, . . . , X„} 
< c(£;|{Xi,...,X„}A{Xi,...,X„}|) 



n 

{d-p)/d 



^-{d-p)/d 
-{d-p)/d 



n 



-{d-p)/d 



(d-p)/d 



< C{df"''K) 



n 



-{d-p)/d 



^^-P'^/'^^-Pid-p)/d_ 
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Therefore, 



J[0,l]d 



Now, we choose m so that {nm ^/'^ — m ^^'^ p^^'^, i.e., we choose m — ?T,'3(d-p)+d. 
Then, Theorem 4 follows from the above bounds for the terms I, II, and III. | 

It is an interesting problem to generalize Theorem 4 to cover more general contin- 
uous probability density functions /. For example, in the case / e [0, 1]*^), 
1 < /3 < 2, we may work with a degree-one block density function; a probability density 
function of the form 

i=l 

where > is a degree-one function and where {Qi, 1 < i < m^} is a partition 

of [0, 1]'' into m'' subboxes of edge length m"^, is a degree-one block density function 
of level m. For this degree-one block density function using Theorem 4 we can get a 
rate of convergence similar to Theorem 3. If we can get an approximation result using 
a degree-one block density function similar to Lemma 7, we can extend Theorem 
4 to cover the case 1 < P < 2 (and hopefully by iterating this argument) and the 
case 2 < /? < oo. However, we face some technical problems in this approximation 
procedure. We leave this problem to the interested reader. 
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